Conditions for non-monotonic vortex interaction in two-band superconductors 
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We describe a semi-analytic approacli to the two-band Ginzburg-Landau theory, which predicts 
the behavior of vortices in two-band superconductors. We show that the character of the short-range 
vortex- vortex interaction is determined by the sign of the normal domain - superconductor interface 
energy, in analogy with the conventional differentiation between type-I and type-II superconductors. 
However, we also show that the long-range interaction is determined by a modified Ginzburg-Landau 
parameter k* , different from the standard k of a bulk superconductor. This opens the possibility for 
non-monotonic vortex-vortex interaction, which is temperature- dependent, and can be further tuned 
by alterations of the material on the microscopic scale. 
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I. INTRODUCTION 

Multi-gap superconductivity arises when the gap am- 
plitudes on different sheets of the Fermi surface are rad- 
ically disparate, e.g. due to different dimensionality of 
the bands for the usual phonon-mediated pairing, as is 
the case in MgB2,^ or due to the repulsive pairing inter- 
action, as it appears to be the case in recently discov- 
ered iron-pnictides<^i^ The other examples of multi-gap 
materials include OsB2, iron silicides such as Lu2Fe3Si5, 
chalcogenides (NbSe2), but also the conventional super- 
conductors such as Pb when reduced to nanoscale4 

In a strong magnetic field all superconducting conden- 
sates form normal-metal voids, as an intermediate state 
before superconductivity is fully destroyed. These nor- 
mal domains tend to merge in type-I superconductors in 
order to minimize their positive surface energy, whereas 
in type-II superconductors they have negative surface 
energy and split into quantized vortices. However, in 
2005 Babaev and Speight predicted the so-called semi- 
Meissner state in two-band superconductors,^ the state 
with localized regions of high and low vortex densities, 
arising from short-range repulsive while long-range at- 
tractive vortex- vortex interaction. This vortex behavior 
was recently visualized by Moshchalkov et al.,^ in the 
form of stripes and clusters of vortices in a single-crystal 
MgB2. Such vortex configurations stemming from the 
long-range attractive vortex behavior (see also Ref. 7 
for review) are clearly very important in the field of su- 
perconductivity, but they also present a bridge between 
solid-state physics and soft condensed matter, where sys- 
tems with competing interactions are of abiding interest.^ 

To date the matter of competing vortex interactions 
in two-band superconductors has not been conclusively 
settled although recent years saw a surge of activities 
in this field. The original prediction in Ref. 5 concerns 
only the case when one band is type-I and the other type- 
II, although it is unclear how different types of behavior 
between bands in k-space (not real space) can be dis- 
cerned. Ref. demonstrated such vortex behavior in 
systems where just one band is fully superconducting. 



and the other superconducts only due to direct coupling. 
Dao et al. found different types of possible vortex- vortex 
interactions and several resulting exciting vortex config- 
urations, but did not provide a universal criterion to a 
priori determine the type of vortex interaction^^ Finally 
some authors expressed scepticism to nonmonotonic vor- 
tex interaction; Geyer et al. showed that the normal 
metal/two-gap superconductor surface energy close to Tc 
depends just on a single Ginzburg-Landau (GL) param- 
eter K, and thus only either repulsive (type-II) or attrac- 
tive (type-I) vortex-vortex interaction is possible. This 
point was later reenforced by Kogan and SchmalianJ^ 



II. METHODS AND DERIVATIONS 

A. The Ginzburg-Landau formalism for two-band 
superconductors 

In this paper we derive criteria for the appearance of 
non-monotonic interaction of vortices in two-gap systems 
described by the standard GL model. Our analysis is 
based on the two-band GL theory, but with correct micro- 
scopic parameters obtained either from theoretical band 
structure calculations or by fitting the experimental pen- 
etration depth or specific heat data by the so-called 7- 
model^^. We begin from the GL energy functional, which 
comprises single-band contributions from both conden- 
sates, the coupling term, and the energy of the magnetic 
field in and around the sample: 
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Here the two Cooper-pair condensates are described by 
the order parameters ^'i and 5*2, H is the applied mag- 
netic field and h the net one. The Josephson coupling 
term provides the 'minimal coupling', well described in 
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literature. The temperature enters the energy expression 
through aj=i.2, hnearly dependent on the temperature 
term r = In T^/T w 1 - T/TcM- The expansion leading 
to the GL theory is strictly valid only in the immedi- 
ate vicinity of T^, but we use this theory at somewhat 
lower temperatures as well, arguing that GL theory qual- 
itatively well describes important physics away from Tc 
(as was demonstrated at many prior instances). Finally, 
it was shown in Ref. lH, that standard two-band GL 
theory contains incomplete terms that estimate ip with 
precision to r^/^. The authors reduce the theory by elim- 
inating latter terms, which results in a single coherence 
length for both order parameters of a two-band super- 
conductor. This is however not a correct physical pic- 
ture at low temperatures, and two coherence lengths for 
the two-band superconductors can be recovered even in 
the GL domain in the extended model of Ref. [l^ Un- 
fortunately, the latter model is presented in the absence 
of magnetic field. To be able to capture all the essen- 
tial physics, at least qualitatively, we base our study on 
the compromise standard GL model. Note however that 
our further explained semi-analytic approach can be ap- 
plied to any improved form of the energy functional for 
two-band superconductors. 

We next calculate the vortex-vortex interaction in a 
similar fashion to Ref. but within a correct microscopic 
framework. The parameters in Eq. ([l} can then be ex- 
pressed as: aj = —N{0)njXj = ~N{0)nj{T — Sj/njrj), 
= N{0)nj/W^, rrij = 3W^/N{0)njv] and F 

iV(0)Ai2/77, where A = , 

matrix with determinant r/; Uj {N{0)) denotes partial 
(total) density of states, vj are the Fermi velocities in 
the two bands, and = 8tt^T^ /7C{3). For details on 
constants Sj we refer to Ref. [ij. This allows us to tech- 
nically define the coherence lengths = -y=^ and pen- 
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parameters Kj — of the two condensates, as if they 

were independent. These are however just parameters of 
the model, and are related only indirectly with the re- 
sulting penetration depth and the healing lengths of the 
two order parameters in the two-band material. Notice 
also that ai and a2 change sign at different tempera- 
tures. In particular, close to Tc both aj are positive but 
the coupled system is still superconducting. Such situ- 
ation is already different from the one studied in Ref. 

where at least one aj was negative. The Ginzburg- 
Landau equations minimize the functional from Eq. ([l}, 
and read (in dimensionless form) 



where jj = 3? [**(-iV - A)*j] , a = {vi/v2f, 7 = 
r/ni7V(0), both order parameters are scaled to W, dis- 
tances to ^1, and vector potential to /ic/4e7r^i. 



B. Long-range vortex interaction 

In what follows, we demonstrate the method to deter- 
mine the asymptotic long-range interaction of vortices, 
before going into fine details at short vortex-vortex dis- 
tances. In cylindrical coordinates, considering the ansatz 
for one circular symmetric vortex ^ j — e^^ fj{r), and sub- 
stituting the gauge A — a{r)9/r, we rewrite Ginzburg- 
Landau Eqs. ([2ti-c) as 
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For r — > 00, a converges to 1 and fj to a constant aj. 
The limit r ^ 00 leads to the set of non-linear coupled 
equations for a^: 



(Xi - ai)ai + 702 = 0, 



a{x2 - 02)02 -I |ai = 0. 



(4a) 
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These can be decoupled by defining the ratio p — ai/a2, 
which then obeys the fourth order equation 
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Such an equation has a laborious analytical solution 
known as Ferrari's method, which will not be presented 
here, but can be found in Ref. [l^. From Eq. (U), one 
obtains the dependence of the constants aj on the ratio 
p as 
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i-iV - Ay^2 - a(x2 - |*2r)*2 - - 0, (2b) 
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In order to eliminate high order terms for large dis- 
tances, we must use auxiliary functions that approach 
zero as r — > 00, namely, Q{r) = o(r) — 1 and (Tj{r) = 
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first order terms in tliese func- 



1 — 3ai) fJi + 7(72 — 0, (7a) 



3a^)a2 + -^ai =0, (7b) 



fj (r) — Uj . Keeping only 
tions, Eqs. ^ become 

(Per I ^ 1 dffi ^ . 
dr'^ r dr 
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wliere we defined A ^ = (ai/Ai)^ + (a2/A2)^. Tlie solu- 
tion of Eq. (I?!;) is the Modified Bessel function Q{r) — 
SsrKiir^i/X). Similarly, if 7 = 0, Eqs. ^) and dZb) 
are decoupled and easily identified as Modified Bessel 
equations, whose solutions are (Ti(r) = riiKo{^/2xir) 
and (72 (r) = ri2Ko{\/2ax2r). On the other hand, if 
7 7^ 0, the equations for aj are still coupled and, in 
order to decouple them, one must define the operator 
L2 = + a (x2 - iaVj, so that L2cr2 = -{jKl/anDai, 
and apply it on Eq. ([7^), obtaining 

V^V^cri + CiVVi + C2(Ti = 0. (8) 

Here Ci — (xi ~ 3af ) + a (x2 — 303) and C2 = 
a (X2 — 3a2) (xi — 3ai) — j'^Kl/anl. The operator 
for axially symmetric solutions has eigenfunctions given 
by Bessel functions Jo{fir) and lodS?"), with eigenvalue 
— or modified Bessel functions Io{l3r) and iiro(/3'"), 
with eigenvalue /3^. From these four eigenfunctions, 
only the latter satisfies the condition that aj must de- 
cay monotonically with r. Substituting V'^KQ{f3r) = 
P'^Ko{Pr) in Eq. ([8|), one obtains 

+ Ci/?2 + C2 - 0, (9) 

and 

cTi (r) = (5i cos{uj)Ko{P-r) - (52 sin(w)ifo(/3+7-), (10a) 

(72(r) = ,5i sin(w)ifo(/3-r) + 82 cos (a;)i^o (/?+?•), (10b) 
where 

Notice that in Eqs. ([TUl) . each must contain the Bessel 
functions for both in a combination that is conve- 
niently written in the form of a mixing angle In the 
7 — >■ limit, one has /3_ -> V^Xi and /3+ — >■ ^2ax2- 
Moreover, substituting Eqs. (fTOl) in the differential equa- 
tion ([7^), one obtains 

^ /3i + (x'!-3a?)' 



so that 7 — > leads to w and, consequently, to 
<Ti(r) ^ 77iifo(V'2xir) and cr2(r) ^ 772^^0 (V2aX2?'), as 
expected. 

The parameters 5k {rjk) in the expressions for Q{r), 
CTi (r) and a2 (r) are unknown real constants that can only 
be determined by fitting numerical solutions for Eqs. ([3]) 
in analogy to what is done in Ref. [l3- 

Having the asymptotic form of the order parameters 
and the vector potential, we now follow the standard 
procedure^^ for finding the vortex-vortex interaction in 
the r — >■ 00 limit, obtaining 



E2B{r) = SIK, Q SlKo - SjKo , 

' (12) 

where the units are now explicitly shown. Here, we list 
the consequences of the above asymptotics. i) Comparing 
Eq. ([T2I) to the one-band case^^, where 



SiB(r) - SlKo (r/Ais) - SIKq (V2r/aB) , (13) 



shows that the lengthscale A~^ = (ai/Ai)^ + (a2/A2)^ is 
playing the role of an effective penetration depth for the 
two-band superconductor in accordance with Eq. (60) in 
Ref. il contrary to A^^ ^ {1/Xi f + {1/X2)^ used in 
Refs. la and [21I which holds only in the (unrealistic) ab- 
sence of coupling, ii) The parameters Sk are in general 
different from each other, but can be calculated exactly 
in the Bogomol'nyi point for the two band system as 
5'f — $2 = 2(5|. For 7 = 0, the choice of ^ = ^ = 1 and 
Ki = K2 = 1 in the two-band case is thus analogous to 
the Bogomol'nyi point kib = l/v^ for the single band 
case and, accordingly, the long-range interaction must 
vanish (and change sign for ki = k.2 < 1)- This directly 
illustrates that coupling of two (nominally) type-H con- 
densates may lead to a type-I behavior of the coupled 
system! iii) In Eq. ([TS]) for single band superconductors, 
it is clear that if kib = Xib/^ib > l/v^ (< the 
interaction potential Eib (r) will be repulsive (attractive). 
For two-band superconductors, Eq. shows that the 

relevant parameters are k1 — -7^-7-, rather than the 

v2 5i 

nominal GL parameters kj for each condensate. If ei- 
ther K,*^_ or K*i are below l/-\/2, the long-range vortex 
interaction is attractive (type-I like). Eqs. Q and (jlOcp 
provide simple means to evaluate this condition, iv) In 
the presence of coupling, the long-range behavior of both 
CTj depends exponentially on the smallest of /3_ and /3+ . 
Therefore, in the coupled case, we can define not only 
a single penetration depth for both bands, but also the 
order parameters for both condensates exhibit the same 
decay at large distances which implies a joint coherence 
length r =ei/min(/3+,/3„). 
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C. Surface energy and the short-range vortex 
interaction 



^2 „,^2 

= - a (X2 - ^-a) ^-2 - 



^ ^2 1' 



(16b) 



The analysis in the previous subsection brings us to the 
discussion of the real criterion for the attractive / repulsive 
nature of the vortex interaction. In the single-band case, 
the changing sign of the normal domain - superconductor 
surface energy Eg at the Bogomol'nyi point is a correct 
criterion. However, in the two-band case and for large 
vortex-vortex distance, the Bogomol'nyi point is deter- 
mined by a single valued k* = min(«;^,Kl) = 
which is not necessarily where the surface energy of the 
normal domain (vortex) changes sign! 

The sign of the energy of the interfaces between 
normal-metal domains and the superconductor deter- 
mines whether merging of those domains is energetically 
favorable or not (i.e. if the superconductor is type-I 
or type-II). In the case of vortices, the smallest possi- 
ble normal domains, the positive vortex-superconductor 
surface energy therefore means that the vortices should 
repel (at least at short distances) in order to avoid the 
formation of a giant vortex. We here show how to cal- 
culate the normal-superconducting interface energy and 
by that predict the type of the short-range vortex- vortex 
interaction. 

We follow a similar approach to that of Ref. [Ill, but 
we take into account the Josephson coupling and the 
temperature dependence of the Ginzburg-Landau (GL) 
parameters, within a correct microscopical framework. 
Namely, we consider the interface between normal and 
superconducting region as the yz-plane at a; = and cal- 
culate the surface energy Es using the one-dimensional 
GL functional at the thermodynamic critical field Hcc, 
which reads 

Es= f dx{2{^f + A^-^l) + {2x1- 

J —OO ^ 

K2 

-2j^i^f2 + (h,, - 72^1 A') ' } , (14) 



A" 



*2 



A. 



(16c) 



The boundary conditions in the normal state {x — > —00) 
and deep in the superconducting state {x — > 00) are 
^]-j{x —00) — 0, A'{x — > —00) = 1, Vj(a; 00) = 
and A'{x 00) = 0. 



D. Constrained GL equations for fixed vortices 

We supplement our argumentation by numerically ob- 
tained vortex-vortex interaction potentials (in a simi- 
lar fashion as in Ref. [22| ). Since the problem of two 
vortices does not have circular symmetry, we now con- 
sider the fixed-vortex ansatz in Cartesian coordinates 
'^j = e™^^^e™^^^ fj{x,y), describing two fixed vortices 
with winding numbers ni and 712, where e™'=^'= is written 
in Cartesian coordinates as 



Xk + iVk 
Xk - iyk 



(17) 



and fk = {xk,yk,0) is the in-plane position vector with 
origin at the center of the vortex k. For the case of two 
vortices separated by a distance d, we take ri = (x — 
d/2,y,0) and r2 — {x + d/2,y,0). With this ansatz, 
the Euler-Lagrange equations for the energy functional 
in Eq. (1) read (see also Ref. 22) 



V2/2 



2 — 2 

X +Y 



-2{A,Y 



AyX) 



A' 



+(xi-./i)/i+7./2 = 0, (18a) 



+ y'^ + 2{A^Y - AyX) 



A^ 



/2 



2 

+a (X2 - /2 ) /2 + ^/i = 0, (18b) 



ant 



where the gauge potential is chosen &s A — (0,A(x),0) 
and ^ j=i.2 are taken real. The thermodynamic critical 
field of the coupled system Hcc is obtained from the con- 
dition that the GL functional in Eq. ([1]) converges to zero 
for H = Hcc, leading to 



Hlc^H\i)al{2xi 



-H^ "2 



c(2) 



£12(2X2 



)+47-f^c(i)«ia2- 
. 

We then find '^i and A that minimize £^5 by numerically 
solving the set of Euler-Lagrange equations for the func- 
tional in Eq. (|14p , which are exactly the one-dimensional 
versions of Eqs. (H^-c): 



j{2 

*i - - (xi - *i - 7*2, 



(16a) 



and 

V X V X i' = 



where 



rl 



A 



ri 



712X2 

'2 



?^2^2 

r2 



(18c) 



n2y2 

^2 



and the angular unit vectors around each vortex are writ- 
ten as 0k = i-yk/rk,Xk/rk,0). 

Eqs. (|18b .-c) are thus the GL equations for the two 
fixed vortices, and we solve them numerically by a relax- 
ation method. The obtained order parameter and vector 
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FIG. 1: (color online) The normal domain - superconductor 
surface energy Es as a function of the ratio of the density of 
states in the two bands (a) and the corresponding vortex- 
vortex interaction energies (b-d) for indicated parameters. 
The short-range interaction force changes sign when the sur- 
face energy changes sign. 



potential are then substituted back in the energy func- 
tional, yielding the energy E{d) for the vortex pair at 
distance d. Repeating this procedure for different vortex- 
vortex separation, we obtain the interaction potential 
AE = E{d) — E{Q) between vortices in the two-gap su- 
perconductor, as shown in Fig. [TJb-d). 



III. RESULTS AND DISCUSSION 

We now apply the techniques described in the previ- 
ous section to calculate (i) the asymptotic long-range GL 
parameter k* , (ii) the normal domain - superconductor 
surface energy Es , and (iii) the full vortex- vortex poten- 
tial (using the constrained GL equations). 

As a first example. Fig. [T] shows the surface energy 
Es and the numerically obtained vortex-vortex interac- 
tion potentials for a set of parameters corresponding (ar- 
guably) to MgBa: ki = 3.71 and ^1/6 = vi/v2 = 0.255 
are taken from Ref. 0, the coupling matrix is obtained 
from Ref. [l^, the temperature is fixed at T = 0.82rc, 
while we vary the density of states in the two bands. We 
note that in all considered cases k* < l/-\/2 and the long- 
range interaction is always attractive, whereas short- 
range interaction changes to repulsive exactly when the 
surface energy Es changes sign with increasing ni/n2. 
To conclude, the long-range vortex-vortex inter- 
action is determined by k* with respect to l/\/2, 
while the short-range behavior is determined by 
the sign of the surface energy Es- This also proves 
insufficient the initial premise in Ref. |6) that if the system 
has A/^i > I/V2 and < l/"\/2, the vortex interac- 
tion should be long-range attractive and short-range re- 
pulsive. The actual behavior is far more complex, and 



can be exactly determined as explained above. 

Recent calculations have shown that as T — Tc, only 
type-I or type-II vortex behavior can be observedj^ In- 
deed, by analyzing k* and the sign of Es at T Tc as 
explained above, we always found the same type of in- 
teraction in either long- or short-range limit. However, 
for T immediately below Tc the sign change of Es and 
the transition from k* > 1/V2 to k* < l/\/2 occur for 
different sets of parameters, opening up the parameter 
space for observation of the non-monotonic vortex inter- 
action. This is shown in Fig. [21 where we plot k* as 
a function of temperature and also the ratios between 
the Fermi velocities (Fig. [5^) and the partial density 
of states of each condensate (Fig. [2b)- The black line 
in Fig. [5| denotes k* = l/\/2 and the white line indi- 
cates where Es = 0. At T = Tc these lines coincide, in 
agreement with Ref. [Ill, but as T decreases, the lines 
separate, bordering the region where the system exhibits 
short-range repulsion {Es < 0) and long-range attrac- 
tion {k* < I/V2), i.e. non-monotonic vortex interaction. 
This finding further creates a new possibility of tuning the 
magnetic interactions in two-band superconductors by 
changing temperature. For example, for the parameters 
of MgB2 given in Ref. ^ (Fig. ^ for vi/v2 = 0.255), we 
find that non-monotonic vortex interactions occur only 
for T < 0.49 Tc— whereas pure type-II behavior is ex- 
pected at higher temperatures. The experiment in Ref. 
6: was done at T « 0.1 Tc, and could thus be repeated at 
higher temperatures to verify our prediction. 

In Fig. [3l[a) a similar phase diagram is constructed 
for recently discovered, and for many reasons exciting, 
pnictides. In particular, we show the results for LiFeAs, 
using the parameters given in Ref. H, except for the 
fact that A12 in the A matrix must be taken negative 
due to the s± pairing. For this material, we extract 
Ki = 2.4, ni/n2 = 1.384 and Vi/v2 = 0.722. Inter- 
estingly enough, as K2 = ki \fri\v\Jn2^i\., we note that 
both nominal GL parameters of the bands are larger 




0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 



FIG. 2: (color online) The long-range interaction phase dia- 
gram for a MgB2 crystal, at different temperatures and for 
varied values of the ratio between (a) the Fermi velocities, 
and (b) the partial density of states, of the two bands. In 
each panel, the black line separates the regions of long-range 
attraction and long-range repulsion (left to right). The white 
lines indicate where Es changes sign and the short-range in- 
teraction changes from attractive to repulsive (left to right). 
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FIG. 3: (color online) (a) The long-range vortex interac- 
tion {vi/v2,ni/n2,T) phase diagram for LiFeAs, for other 
parameters taken from Ref. Si. The shown isosurface corre- 
sponds to «:* = 1/ \/2 and the change of the long-range vortex- 
vortex interaction, (b-d) 2D cuts of (a) in the T = 0.9 Tc, 
vi/v2 ~ 0.722, and nxjn^ — 1.384 planes, respectively. Black 
(white) lines correspond to k* = l/\/2 {Es = 0). 

than l/-\/2 if y/n^vf/n2V^ > 0.295. Therefore, it can be 
once more verified that in a large portion of the parame- 
ter space where both bands are convincingly type-II, the 
coupled system exhibits type-I behavior. In Fig. [3][b), 
we show that for T ~ 0.9 Tc, the Es — (white) and 
* = l/y/2 (black) curves coincide for small ni/n2 and 
large Vi/v2- This behavior persists even at lower temper- 
atures, as shown in Fig. [3l^c). However, in the opposite 
case (large ni/n2 and small V1/V2), the curves separate, 
forming a region of non-monotonic vortex interaction in 
the phase diagram which grows larger as temperature 



decreases (see Fig. |3l^d)). This broad temperature range 
for the observation of partial vortex attraction is impor- 
tant experimentally, to discriminate the non-monotonic 
vortex interactions from irregular vortex lattices formed 
due to intrinsic defects in the materiaP^ (with latter be- 
ing dominant only at temperatures where the vortex core 
and the defects are similar in size, unless defects are of 
magnetic nature). 



IV. CONCLUSIONS 

In conclusion, we have demonstrated the semi-analytic 
method to relatively easily determine the nature of 
vortex-vortex interaction in two-band superconductors. 
This is of significant theoretical and experimental impor- 
tance, as Figs, m 13] sketch just two examples of many 
possibilities attainable by two-band hybridization. Note 
that a plethora of transitions, even reentrant behaviors, 
can be found as a function of the microscopic param- 
eters, which can be tuned experimentally (to some ex- 
tent) by e.g. carrier injection. ^5, Finally, with appro- 
priate modifications of the initial energy functional our 
approach can also provide insight in similar situations 
encountered in nanoscale superconducting films, tailor- 
made two-component superconducting hybrids, and dirty 
two-band compounds. 
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